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Biographical data will not be published unless there 
is a demand for it. This is due to pressure for the pub- 
lication of material which seems to be of greater interest 


to our readers. 
Editor. 


EDITORIAL NOTE - 


Readers are invited to send us their selections of 
remarkable situations in Mathematics that can be explained 
at the freshman level. 

H.V. Craig, Editor. 
(Answers, continued from page 57) 
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THE RANDOM SIEVE 
by David Hawkins 


Cramer (1) discusses a simple stochastic model for the distribution 


of primes. Let there be a sequence of independent trials of an event 
with 


(1) Pr(S,) = 1/log n 


The numbers P,, P,, «.- Pl, «++, of which P is the mth value of n 
for which Ss, occurs, will then with probability one have a limiting 
density 1/log n, like the primes. A number of other consclusions follow 
from the strong law of large numbers, for example that with proba- 
bility one 


P 


P 
(2) li 


The model is however quite artificial in that it has the prime 
number theorem built into it ad hoc and in its assumption of inde- 
pendence. Primes are not independent in a statistical sense. The 
occurrence of an unusually long run of composite numbers leads one 
to expect a compensating increase in the number of primes later on. 

A more natural stochastic model is that of the random sieve. In 
the sieve of Eratosthenes we sieve out the multiples of every number 
which is not,a multiple of some earlier sieving number. We define 
the random sieve as follows: Check the number 2, and then with pro- 
bability % strike out each subsequent number. If P, is the first 
number not stricken out, check it and strike out each number there- 
after with probability 1/P,. P, is the next number not stricken out, 
and we use 1/P, as the probability with which to strike out each 
subsequent number, etc. The set of all possible sequences of numbers 
checked (sequences of random sieving numbers) contains the sequence 
of primes as one of its most probable members. Another typical member 
is the sequence of ‘‘lucky numbers’’ (2). S_ now stands for the pro- 
position that n is a sieving number, and we have the following 
recurrence relation: 


Pr(s,)? 


n 


(3) Pr(S, , ,)=Pri(S,) - 


Proof: Let T, be the contradictory of S . Let S% be the proposition 
that n is not sieved out by any sieving number less than n - 1, and 
let T* be its contradictory. Then 
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Pr(T.T,,3) = Pr(T,)?, Pr(S,T%,,) = Pr(S, Pr(T,), 


But 


and 
Pr = Pr(S,)?/n. 


Since the alternatives are mutually antasie, the recurrence relation 
follows. 


We solve (3) as follows: 
Substitute g, = 1/Pr(S_), obtaining 


(4) 
Bn +1 1/g, 


which may be solved recursively, using g* = 5 1/n =L(n) as the 
first trial solution. This gives g, = L(n) + 0(1) and thus 


(5) Pr(S,} ~ 1/log n 


Thus the random sieve gives asymptotically the result assumed 
ad hoc by Cramer. The events S are not, moreover, independent. In 
fact, it is obvious that 

(46) Pr(T T. 1 eee 


n+ 


p-1) = Pr(T,, yr 
which shows that the interdependence is negative, in conformity 
with our comment about the primes. The result does not affect the 
validity of (2), however, unless it strengthens it. For if in (4) 
we put r= e(leg 2n), it is easy to see that the occurrence of a 
run of that length without sieving number is asymptotically of 
probability 1/n°. From the convergence of 51/n° for c>1 it follows 
that with probabilicy one the number of runs of such length is 
finite. Hence log 7, is almost certainly an upper bound, from 
some on, to the interval Pq: Because of interdependence 
it is more difficult to prove that this is false for any c > l. 

It is not difficult to define and solve various random-sieve 
problems analogous to those of multiplicative number theory, for 
example the relative frequency of pairs of sieving numbers separated 
by a given interval, the expected number of “divisors” of any 
number, etc. 

In the case of the lucky numbers 1, 3, 7, 8, 13, 15, ..., we 
start with odd numbers and sieve out first every third odd number, 
leaving 7 as the next lucky. Then we sieve out every seventh of 
the remaining numbers, etc. In this case nothing is known about the 
asymptotic density, so we are in the position, say, of Gauss viz 
a viz the prime distribution. A randomization of this process is, 
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however, just the random sieve again, and from this fact we im- 
mediately conjecture that the asymptotic density of luckies is 
1/logn rather than, say, B/logn where B 7 1.* 

There is one direction in which the random sieve may facilitate 
something more than conjectures. The distribution of numbers prime 
to the first m primes, or that of pairs of such numbers separated 
by a constant interval (e.g. 2) is in all likelihood more regular 
in a certain sense than the corresponding random sieve distribution, 
and if this is true it implies a number of results somewhat stronger 
than those that have been obtained, e.g. the infinity of twin 
primes. 

Suppose that the first m random sieving numbers are the first 
m primes, PL = =p, r=1, 2, ... m. Then the distribution of the 
number of numbers not sieved out from a sequence of NV consecutive 
numbers > p, is given by the binomial distribution, with proba- 
bility Q = , = 1/2)(1 - 1/3) ... (1 - 1/p,). This same distri- 
bution may be expressed in a different way as the sum of 2" random 
variables (non-independent). Let N(i) be the number of numbers 
sieved out by the sieving number p;, N(i,j) those sieved out both 
by p; and Pj» etc. Those not sieved out will be given by the well- 
known combinatoric formula 


1, 


For the sieve of Eratosthenes, on the other hand, we have a 
precisely similar expression. In this case the distribution of 
numbers prime to 2.3.5. . p, =K, is periodic modulo K,, and we 
pick our sequence of N numbers, N << K,, at random from a period 


of length K,. We can calculate moments for the distribution of the 


number given by (7), both for the random and the Eratosthenes case. 
If as seems likely we can prove that moments of even order for the 
Eratosthenes distribution are smaller than the corresponding moments 
of the random sieve distribution, then it will follow that the 
longest interval between numbers prime to K, is of the order of 
Pp, log p,, a stronger result than has been obtained by other 
methods. A similar argument applied to twins prime to K, would, if 
valid, sstablish the infinity of twin primes. It is net hard to 
prove the inequality for second moments, but the problem of proof 
for moments of order 2k remains. 
*e @ 

* Subsequently verified by W.E. Briggs and the author, and inde- 
pendent ly by P. Erdos. See *The Lucky Number Theorem" to be published 
in this magazine. 
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A Slide-Bee 
James G. Dyhikowski 


Recently there appeared in this magazine an article entitled 
“Dig That Math:’ The article was humorous and somewhat sarcastic. 
It implied that mathematics would need a radical change before 
mathematics teams would he formed. I would like to quote an article 
date lined Chicago, March 19, 1957 and entitled,“Mathematics is a 
New ‘Sport’?’ 

“A group of engineers who are worried about the nation’s shortage 
of engineers have devisedaslide rule competition for youngsters. 
They hope it will ‘make mathematics as popular as sports.’ 

“The first of a series of planned ‘slide-bees’ was held yesterday 
in suburban Wheaton High School. Seven-member teams from seven high 
competed, each armed with a slide rule. 

“The contest was run off like a spelling bee, with contestants 
given mathematical problems and those answering them fastest piling 
up points towards the championship. 

“Among questions posed was this one: 


How soon will the Russians overtake our supply of en- 
gineers, if they turn out 81,000 engineers every y ear 
to our 28,000 engineers? Assume that we now have 542,000 
engineers, against 396,000 engineers in maauedl and that 
deaths and retirements are disregarded. 


‘* The answer to this was 4 years and 7 months. 
“The winner was George Guerin of suburban Hinsdale High School. 
**Guerin’s prize: A slide rule. 


The only ambition which I have, is that this movement well spread 
over the whole country and get students interested in science and 
engineering. 
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CALCULATION OF A COMPLETE SYSTEM OF TENSORS 
WITH THE AID OF SYMBOLIC MULTIPLICATION 


Lewis Bayard Robinson 


Introduction 


A tensor is a generalized covariant. Camplete systems of co- 
variants are given by the solutions of complete systems of linear 
partial differential equations. 

Complete systems of tensors are given by the solutions of dif- 
ferential equations of the Riquier-Saltykov type which contain 
complete systems of linear partial differential equations as a special 
case.°* 


Most of us know that the German geometers write a binary form 
thus: 


(ax)? = (ax, + +a,x,) = a, 2a, + 


We shall use the above symbolic multiplication to help usin- 
tegrate the differential equations which arise. 

Dr. T.C. Doyle has written an important paper on the tensors of 
projective diflerential geometry.3- 

The author is extending the work of Wilczynski and his pupil 
E.B. Stouffer the covariants of linear differential equations to 
tensors and like them he rests on the work of Lie but he uses the 
German symbolic notation as an auxiliary. 

Professor Cramlet has written a work: The Invariants of an N - 
ary Q-ic Differential Form.4* But his definition of a complete 
system of tensors differs altogether from that given by the author. 
Nor does he follow the method of Lie and he does not refer to Riquier. 
He bases his work on Aronhold. It is only remotely related to the 
work of the author. 

The first paper published by the author on the above subject 
appeared in the Doklady of the Russian academy (1937). Professor 
Hlavaty of Praha wrote a review of this note and remarked that the 
author had not used a tensor notation. I have with considerable 
profit introduced a symbolic notation in this work. 


1. See Riquier, Les Systenes D’Equations aux Derivees Partielles, page 502 
2. See Grace and Young, Algebra of Inver iants. 


3. Doyle, Tensor Theory of Invariants for the Projective Differential 


Geoaetry of a Curved Surface. Transactions of the American Mathematical 
Society, 1944. f 


4. See, Annals of Mathematics, 1930 page 134 
5 


| 
| 
| 


MATHEMATICS MAGAZINE (Sept. -Oct. 


Part (I) 


In a previous work the author has computed a complete system of 
semitensors associated with the system of diflerential equations 


(s) 


+ + Po2¥d + + = 


In the following paper he will compute a complete system of ten- 
sors associated with the above system. 

The tensors we shall consider are of order r = 1 when there is 
aly one subscript. Such a tensor is I;. I; is a tensor of order r=2. 

We shall now set up the system of equations defining the tensors 
when r = 2. It will then be easy for anyone to write the system 
when r = 1. 

When we were computing a complete system of semitensors the de- 
terminants of the transformation was 


a 


Ory 


12 | 
22 


The determinant of the transformation of our tensors is written 


2 0 
0 
The infinitesimal transformations of A are given by 
1+¢,,5t dt 0 
0 0 1+é'8 


The transformations of the tensors are written when r=1, thus: 


3 
T,= (i = 1,2,3) 


When r = 2, they are written 2 


1. See Robinson, Un Systeme De Riquier et Le Calcul Tensoriel 
Bulletin de la Societe Mathematique, 1940, page 129 


2. The symbol @ indicates symbolic multiplication and I; 


6 
| | 
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+A, AS; + A, + A, (434133 


We shall compute our tensors by Lie’s method of infinitesimal 
transformations. The infinitesimal transformations are: 


+ 21, - + 20 


ST ~ pdt + Ty + - Ty * 


* 21,50, - 21 21596 St 


And from the above follow the formulae denoted by (F).! 


+ (1, "St = — 1,,6,,5t + + 


1, + 21,4 ‘dt et cetera. 

We write the equations analogous to system (A) in the anterior work: 
(I,,) = 0 

(i,j,k, t=1,2) 


= - ®,,(1,,) = - 122 
0 = + 
+ = + Iy2 


1. Wilczynski uses infinitesimal transformations to c te complete system 
of covariants in his book Projective Differential 


5143 + T + p,,)dt - + 1, 

0 = - Ty, 
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= 


= 


- 2I,, = 
= +2I,, = + 
= 0, = 0; 
= + 2123, ,(I,,) = + 2I,, 
= +123, = + 2133; 
%2,(1,,) = 9; 
= 0. 


We derive the above from formulae (F) by equating the coef- 
ficients of the and 


We know how to write the expressions, Q; ,(F), ,(F), 
for the author has given them in the anterior work. But we have 
introduced a new symbol o,(F), which we write 

oF 
5, P12 


where @ = e, 


Wilczynski has computed this seminvariant. See Projective Dif- 
ferential Geometry, page 110. 

We observe that of the above equation Wilczynski uses only 
the terms 


since he does net for his purposes require the other terms, when 


(Sept.-Oct. 
= +2 + | 
= + = | 
= 0 2145 
| 
| 
20 
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‘he is computing only invariants. ! 
And to the equations of the above system we add one more set: 


1 ey, 2 2p, » 206 


All the equations of the above system we will denote by (A). 
Since this system (A) represents the infinitesimal transformations 
of the I,, the totality of all the solutions of (A) gives us the 
complete system of tensors. 
Now let 


= 274 + + YPr2 
Z, = + + Y2P22 
Substituting these new variables in ,-:(F) and O,-(F) we obtain 


(F) 


These are the equations which correspond to system {4) in che 
anterior work. Let us now set up the equations which correspond Lo 
system (B) in the anterior work. We write: 


(F) = 0 


o,:(F) + 27 2I + 
Fy, 


oF OF 
®,,(F) - - 


oF OF 
®,,(F) + + a, 


1. See Projective Differential Geometry, page 121, formula (103). We 
replace ©, by © 
2. See the just cited formula. (103). 


| 9 
| 
- 50 
| 72 20 
| 
oF oF 
| oF 
| 
oF oF oF 
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oF 


— I 
12 a> 


+ 27 


aF 
®,,(F) + 2T +I 


13 23 


OF OF ¥ 


—— - - I,, —— = 0 


-1 


Write the new variables 


After substituting these new variables in the above equations 
we obtain 


- Y2 2 x, 12 


(r)} 


I 


22 


. (I) indicates all five equations. 


—-0 
Y,2#2Z,+ 
Y,=Z,+ ja 
OF oF - 0 
oF 
oF 
— — 
OF - 
1 
oF oF 
33 
33 OF + + 21 33 
72%, 2 ay, 
OF 
22 
oF OF I + 
oY, 
OF oF 
+ 13 3 
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From the above five equations we have omitted, the differential 
OF ¥ 
coefficients » because we do not need them to 
Py 
compute tensors. (I) is a reduced form of (B). 
Also ¥, and Y, satisfy 


0 
0 
= 0 


which therefore vanish identically and appear no more in this work. 


If we can solve the five key equations (I) we can compute the 
complete system of tensors. 


Part (II). 


We shall begin by writing the five key equations for the case 
where r= 1. 


When we have solved these we can solve (II) by symbolic multi- 
plication. 

Lastly on page 14 we will write the formula which gives the 
complete system of tensors. 

The five key equations are 


oF 
ay, 


1. (II) denotes the five equations. 


OF 
3 oF OF OF 
oY, T, 3°, 
(IT) 
3 OF OF 
dy, 3 
oF OF or 
OF or OF OF 
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The fundamental solutions of (JI) are written 


a 
=e {y, I, - } 


=e 
=I? 


as we have already stated. (See pages 8 and 9. 
Now 


(41, - Ty T,+ 2 yy +f 

The above is unaltered by all transformations of 
I,° I, 2° 

But referring to foot note page 5, we replace I ;* I, by Ti; 

and write 
- = 6492 - 12 + Ty) = 


We perceive the above is unaltered by the transformations of 


so it is a solution of (I). 
In a similar way we get other solutions of (I). 


In short the transformations of I,* I I,° I, I, *I, and 
Ty Iy2 are tsomorphic. 


To verify the above statements we need only write 
Ae 
Tye Ty = + 2 44,4) + T, + + 


2 
+434, °7, 


- A2 2 2 
+ + 241451743 + + 245 25 +43 155 


etc. 


I 


12 (Sept. -Oct. 
3 3= DT, 
Ty 
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It is evident that these two pairs of transformations are 

isomorphic. It follows that symbolic multiplication using the no- 

tation of Clebsch Gordan gives us the fundamental solutions of (I). 
For as these authors write 


we in our turn will write 
= 
In F;* replace I I; by ;- Then write 
F; °F, = Fi; 


Clearly the F;; are the solutions of (I). (Seetheorem page 34.) 
So we write 


Fy = ( = + 
F,- F, = % - 


2 2 


et{y, I, - y,1,}+ e*D I, = ,, - 


F,+F, = {¥,1, - Y,I,}+ I, {Y,1,, - Y,1,,) 


F, F, =F,, =e DI, DI,=@ D°I 


We have now obtained all the fundamental solutions of (J) by 
symbolic multiplication. The author has proved that these solution: 
are correct by direct calculation. 

We are now in a position to give the result at which we hae: 
been aiming. . 

The integration of system (I) and (II) is equivalent to the in- 
tegration of two systems of Riquier.! Consequently, with the aid 
of our solutions of (J) we easily compute a complete system of 
tensors for the case r = 2. 

The general solution of (B) is written 


OF Faye Fay C,) 


1. See Requier, Les S stemes D’Equations Aux Derivees Partielles, 
Chapter XIII page 502. 
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where ® is an arbitrary function and the C, are the fun- 
damental covariants of system (S) which Wilczynski has 
already calculated. But the system (A) given above represents the 
first infinitesimal transformations éT,. of the tensors Ti 
(See page 9%). Hence our tensors are the solutions of (A).To solve 
(4) we proceed thus.! 
let the %,. be six independent functions of Fy and C,. 

Write 


,, =0 


33 


22 ~ 223 = 


Solve with respect to Fy,, + We get 
rf t2 33 


F;; =¥,,(C,) (i,j = 1.2.3) ji 


where the ¥ are arbitrary functions. Write 


DF, eeeee F,.) 
Tyo I,,) 


Me 


The F;,; are linear with respect to the I therefore ™ does not 
depend on the I. 
Now we solve (E) with respect to the I; ; and get 


j>i L>k 


are the minors of M 


where 
The abae formulae give us a functionally complete system of 
tensors associated with system (S) when r= 2. 
In a similar fashion one can compute a complete system of tensors 
when r= 1. 
1. The process is described by Riquier in Chapter XIII of his book. 


Baltimore, Maryland. 


= 
1 
wal 
M z tj 
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SOME OPERATIONAL METHODS IN THE 
CALCULUS OF FINITE DIFFERENCES 


Joseph Talacko 


1. Introduction. 


Many formulas of the calculus of finite differences are very 
complicated. They can be simplified somewhat by introducing opera- 
tional symbols for numerical interpolation, numerical differentiation 
and integration, for summation and solutions of difference equations. 
Special forms may be derived readily by the use of operators in 
some operational schemes. Whenever we have equidistant intervals, the 
application of these operators enables one to memorize complicated 
formulas in several basic relations, especially for central or mean 
differences. 

In this paper a general operator, which possesses the distri- 
butive, commutative and associative properties is introduced as a 
generalization of calculus of symbols [1], (3), [8]. The standard 
operators of calculus of finite differences are special cases. 
Furthermore, this general operator may be expanded in infinite series 
and the inverse operator defined. The introduction of hyperbolic 
functions and their polynomial expansions are also demonstrated in 
this paper. 


2. The Calculus of Symbols. 


Let f(x) be a real, single-valued, continuous function in a 
closed interval a < x < 6, possessing a continuous differential 
coefficient of order k. Let hy h,, hy, hi; Hy, Hy, be 
real constants . 


Definition: An operator H is defined by the relation 
(2.1) Hf(z) =H, f(x + +h,) +... +H fle +h). 


From (2.1) we get some particular cases: 


If h, = 0, 4; = 0, t=2, 3, ..., m3; we get 


a= 8,. 


A special operator H = 1, does not change the given function and a 
zero operator is defined to be H = 0. Two H operators are equal 
if, and only if, their application to any function gives identical 
results. Let 


15 


MATHEMATICS MAGAZINE 


K f(x) = + ky) + ... +K fla +k), 
(2.2) 
Lf(x)=L, f(x + 1,) + ... + L, ft + l). 
The necessary and sufficient condition for K = L is that 


This operator possesses the distributive, commutative and associa- 
tive properties. It permits the formation of polynomials and com- 
bination of these polynomials. From the definition (2.1) we get 


(2.3) f(x) + glx)} = Hf (x) + 
Furthermore: 
(2.4) (K + L)f(x) = Kf(x) + Lf(x) 


{ K{Lf(x)] = [KL] f(x), 


L(Kf(x)] = (LK) f(x), 


and simply 
KL = LK. 
By induction, we may write, 
H(KL) = (HK)L 
(H + KML HL + KL 


(H+ K)L = + K). 


From the definition of the H-operator (2.1), we may write also an 
equation 


(2.6) H+X =K 


with an unknown operation X, 
Let us define the operation of division by 


(2.7) Hf(x) = glx). 


This operator H ~! possesses the distributive properties, but it is 
not always commutative. 
Thus successive operations are indicated in reverse order. That is, 


(2.8) = glx). 


The operator HH~! gives the unique result 


--Oct. | 
| 
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(2.9) HH7*=1, 


The H-operator may be generalized by introduction of the operator 
H., for which the limit is 


(2.10) lim H, f(x) = lim g(x,r) = g(x) 


and simply 
H = limH,. 


Similarly, we define 


(2.11) (HH) = lim + + +... + 
r 
lim b 
the symbolical infinite series. It converges, if the H-operator is 
the zero operator. 


3. Special Operators in the Calculus of Finite Differences. 
Let tabulated values of the function f(x) be 
flxo), fle), flx,), «.. 
and let the intervals be equidistant: 
= ° h. 
For the operator H, we have 
(3.1) H f(x) = f(x + h) 


which is equal to the well known operator £. For every integer n, 
define an exponential operator 


H"f(x) = f(x + nh). 


For every integer am 


Kfla)=fle + — 
m 


We have the relations: 
K® =H", 


Of course, instead of a rational number n/m, we may have any real 
number » and the sequence 


} 
17 
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by eee, eee 
if we write 
limp. = 


we may write 


lim F(x) = lin + ph) = +uh). 


r-@ 


Because f(w) is continuous in the neighborhood of w = x + uh, we 
may write simply 


lim = 
Special operators, with respect to the interval A are: 
a. The displacement operator E : 
(3.3) E f(x)= flx +h). 
b. The forward difference operator 4: 
A= 
A flx)=flx +h) - f(x). 
c. The backward operator VY : 
.1 
where 
Vf(x) = flx +h) - = -Vflx - 


d. The central-diference operator 5: 


with 


= f( ( 
- 


e. The mean-difference operator 


-% (gh /2, prh/2) 
h h 
uf(x) =% {flx + > + f(x - 


f. The differential operator D: 


D=h-*lim 


a-0 a 


18 
(3.8) 
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This operator gives the differential of the function f(x) with in- 
increment h. 


Df(z)=h =h *lim 
a 


Recalling the definitions of the displacement and other diterence 
operators, we may write recurrence formulas relating differences and 
diflerentials and find some simple formlas of interpolation, numerical 
differentiation, integration and summation. 

We may get other usable operators, by the introduction of hyper- 
bolic function. Let us first introduce the familiar operational 
identity. 

(3.9) E=1+,jaz=e? 


and the fundamental identities involving central differences, hyperbolic 
sine and hyperbolic cosine: 


= ph/2_ prh/2 _ 


§ = 2 sinh (0/2) 


(3.10) 


(3.11) 
uw = cosh (D/2). 
Solving for the operator D in (3.9) we get 


(3.12) 


This operator D may be written in the form 


(3.13) D=2 sinh” *(;) =2In +¥1 +8 7/4), 


or after substitution 
(3.14) D =2 cosh~4u =“1 + 87/4 . 


Many other relations may be found and various algebraic artificies 
for manipulating the basic identities may be discussed. We will show 
the usefulness of some operators which are similar to those of 


Lehmer [4], [5], and Michel [6]. 


p= WEN? + = + @-D/2) 
@ Ar 
D = In(1 +O) >> 
r 
rel 
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4. Some Formulas for Numerical Differentiation 


First, it is possible, taking powers of (3.12) step by step to 
get the well known formula of numerical differentiation of tabulated 
functions f(x), based on (3.9) in the form: 


= rin 


r=1 


The central diflerence formula offers more elegant results from (3.13) 
for n =2m 


p [2 sinh~*(5/2) 
(4.2) 


= §78[] (8/2)? + 
2 3 2 


As a particular case we have for m = 2. 
1 


7 
p*= —8 2 —— 8* ] . 
4 240 


Ifn =2m+1, we can write 
= [2 28+! 
(4.3) 2 3 ® 
Consider the differentiation with respect to 5 of the equation (4.2) 


using the order as power exponent 


d 
(4.4) — = 2m [2sinh™? (8/2)) (1 4 82/4) = 1 
s 


hence 


pr} 
2m 


In particular, we get for example 
D? = 1/4 85 + 7/120 8’ - ...), 


5. Some Interpolation Formulas. 


Using (3.1) we may get the Newton forward interpolation formula, 


20 | 
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for an interval of length h: 


n 
1 >> 
rs 


and similarly the Newton backward formula becomes 


(5.2) = (1 =1 > 
rel 
The Stirling interpolation formula, using relations (3.6), (3.10) 
and (3.11) may be written for n =2a 


(5.3) = cosh 2m + sinh 2a 


The expansion of sinh aD and cosh aD is simple in terms of sinh D. 
Let 
sinh D=t sothat and 


sinh aD = sinh (m* sinh‘ t) = f(t). 
Differentiating twice with respect to t, if f(0) = 0, and f'(0) =a 
we get the diflerential equation 
(5.4) (1 + t7)f"(t) + tf’ (t) - m@f(t) = 0 
Using the Leibnitz theorem, diflerentiating r-times, for t = 0, we 
get a recurrence relation 
Equations (5.4) and (5.5), with initial values of f(0), f’(0) enable 


us to obtain the successive coefficients in the Maclaurin series of 
sinh (aD). This gives us the Stirling formla: 


(m? - m(m?- 1)(m? - 32) 
D+ 


D+ — s sinh®D + ... 


Expressing sinh D in terms of half arguments, we obtain 


2m (4m?- 4) 
5! 


E** = 2m sinh (2/2) cosh?/2) + sinh? (2/2) cosh(?/2) + 
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2m (4m?— 4)(4m 4 ) 


sinh?(D/2) coshW/2) +... + 


4m? 
(D/2) + ... 


finally 
(5.6) 
1)(m2- 2?) 5 m?(m?2 — 1)(m? 2?) 6 
51 wd” + + 


If we let g(t) = cosh(mD) = cosh(m sinh~'t) we get Ressel’s inter- 
polation formula from (5.4): 
f(0), f' (0) = 0 
2 m2(m2-22) m2(m2-2?) (a? -4) 


*2 4! 6! 


or 


D D 
~ cosh (2m- 1) + sinh (2m - 1) 


= cosh (D/2) + “a sinh?(D/2)cosh (D/2) 


‘ 2m(2m -2)(2m - 4)(2m + 2) 


a sinh*(D/2) cosh (D/2)+... 


(2m —1)2m( 2m -2) 
3! 


(2m -1)(2m)(2m -2)(2m - 4) (2m+ 2) 
5! 


sinh?(D/2) 


+ (2m - 1)sinh (D/2) + 


sinh? (D/2) + ... 


Finally 
(m 


in(m - 1)(m -%) 3 m(m -1)(m @1)(m - 2) 
834 
3! 4! 


54 


+ 


m(m -1)(m41)(m-2)(m-%) 
+ 
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6. Inverse Operators A) yw 


The meaning of the inverse operator, as the operator of inte- 
gration is best demonstrated by the equation 


(6.1) f(x) = D7 '¢(x) g(x) dx +-C, 


which is the inverse operation of differentiation 
Df (x) = glx) = hf'(x). 


The indefinite integral D~! leaves an arbitrary constant of + 
gration. So, as we find the operational expressions for nur 
differentiation, we may find operational formulas for nume.. 
integrations and summations. If’ the operators are taken between 
definite limits, the constant of integration vanishes, and we may 
write 


b 
= glx) dx. 
a 


To find Newton’s integration formula, we use the following expansion 
of the operator D; as a new operator (e? ~ 1)7}: 


= 


where as known Bernouli numbers [4], [5]. The operator 
setisfies the equation 


B B B 
A + D a? ale” 1)(D + 


It is a reciprocal operator to 4, It is not the most general operator, 
but all results of the operation (6.2) will be a subset of results 
of the operation Aw! We may write 


B B 
(6.3) 2p - p?+ 
2! 4! 6! 
or 
D-*f(x)=1/h f flx) dx 
B, 


B B 


which is the known Euler-Maclaurin formla. In general we get 
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and finally, the Newton integration formula 


(6.5) 1 1 1 
f(x) dx= + 5 - 


1 19 
+ A?f lx) 790 A? f(z) + 


Similarly, the operator gn: represents the inverse operator, in 


the sense that it nullifies the other operation: 
(6.6) wd fle) =4f lx) + fle - 1) + fe -2)+ ...4+C 
and for definite limits 


- ¥f(b)+f(b-1)+ ... + fla + 1) +%f la). 


The summation formulas may be developed by applying the procedure 
of differentiating with respect to an operator, finding the new 
operator, say y(x), and taking the limits a,b. 


7. Other Extensions of Calculus of Symbolic Operators. 


By the seme method as we have empolyed in sections 2,3 and 4, it 
is possible to define interpolation formulas in one exponential 
operator, by a method similar to the way Lehmer [6] defined the 
Bernoulli numbers. For example, the discussed Bessel’s formula may 
be defined as a primitive function, using hyperbolic function 


cosh(m - %)D 


inh(m - %)D. 
cosh(D/2) 


If we consider the other operator, say (E" -C)", [6], [8], we 
may develop a simple procedure for solutions of diflerence equations. 
If we consider the infinite series in some operators, through func- 
tional expressions, we get a reliable tool with may applications. 
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MISCELLANEOUS NOTES 


Edited by 
Charles K. Robbins 


Articles intended for this department should be sent to Charles K. 
Robbins, Department of Mathematics, Purdue University, Lafayette, Ind. 


A FINITE SEQUENCE AND A CARD TRICK 
Ali R. Amir-Moez 


Suppose we take a number of cards, or we might choose, for ex- 
ample, C blank cards and write a number on each of them with 
repitition allowed. We pick a number a, larger than all of the 
numbers written on the cards. Holding the cards face up, we count 
from whatever number is on the top card, up to a. In this way we 
obtain a pile of cards. We put this pile down in such a way that 
the card from which we started counting is on the top of the pile, 
but not necessarily face up. In fact it is better to have it face 
down. We continue this counting until either we exhaust the original 
pack of cards, or until we get a remainder r such that r <a. With- 
out turning up the top cards we can find N, the sum of the numbers 
on those cards, knowing only C, the number of the cards, p, the 
number of the piles, r, the remainder, and a. 

Let nM, 5 sees ny be the numbers written on the top cards of 
the first, second, ..., p th. pile respectively, and N,, N,, ..., 
N. be the numbers of cards in the corresponding piles. It is easily 

rved that 
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=a-n,+1, 


! 


No=a-n, +1, 


=z 


+8, 
Adding these equalities we get 

N, +N, +... +N, =pa- (mn, +n, +... + + p. 
But my and N,+N,+... +N, =C-r. 


Therefore N=pla+1l)+r-C. 


Now we have the formula. Let us entertain our friends. Pick up 
a deck of 52 cards and let Jexll, Q<=12, and Kz=13. Suppose a=15. 

Ask a friend to do the counting in the way described above in 
the general case. As an example, he may see 3 on the top and he 
counts 3, 4, ..., 15, paying no attention to the numbers on the 
cards after the 3. Let the next card on the top be 6. Then the 
counting will be 6, 7, ..., 15. Next he sees 2 on the top and he 
counts 2, 3, ..., 15. Let the next card be J = 11. He counts 11, 12, 
«++, 15. Now a King comes on the top and he says 13, 14, 15. Finally 
the card on the top is 7, he counts 7, 8, ..., 13 and there are no 
more cards. Of course you are not watching him count the cards, in 
fact, you may be in a different room. He only tells you that he has 
5 piles and 7 cards are left over. Substituting in the formula for 
N we have 


N = 5(16) + 5 - 52 = 35. 


So you say that the sum of the top cards is 35. Clearly the top 
cards in this example are 3, 6, 2, 11, 13, whose sum is 35. 

There is a short cut for getting N. We leave it to the reader 
to find it. 


Queens College, New York. 


(Miscellaneous Notes continued on page 41) 


MULTIPLE NUMBERS 
John A. Tierney and John Tyler 


A complex number is a number of the form x + iy where the com- 
ponents x and y are real and i = V-1. If we replace x + iy by x + ry 
where r? = -1 we may regard x + ry as representing two complex numbers 
x + ly and its conjugate x - iy. It is interesting to consider the 
generalization to n components where r satisfies an equation of degree 
n. 


We define a multiple number as a number having the form c = 
n 


x,r' where the x; are real and r satisfies the equation paris 


n 
a,r', a, real, with the latter equation having no mltiple roots. 
iz0 
c is a multiple number in that it represents simultaneously n + 1 


complex numbers. 
We consider the case n = 1 and take c = x + ry where r satisfies 


r? = pr +q having roots 


2 
p +4q p 
(1) 1’ 2 7 0.) 


We call x and y the components of c and x, + ry, and x, + Try, are 

defined to be equal if x; + = + and + = + 

If x, + ry; =%2 + ry, their components are equal, otherwise the 

equation (x, - x,) + (y, - y,)r = 0 would have two distinct roots. 
The norm N of c is the product of x + rsy and its conjugate 

x or Ne) = x? 4 (r, + + =x? + pry - 

The conjugate of x + ry is x + (p - r)y since their product is N. 


To obtain a polar form for c, let c =x + ry=¥e"? 


6 


y= VN 8 


=x +¢ry =WN e("P/2)8 
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We call YN the modulus of c and @ the amplitude of c. The mul- 
tiple number e{"*?/2)® has unit modulus and we write 


(2) e(r-p/2)® 4 rs(8). 


To find the functions c and s we substitute r, andr, in (2) and 
solve for ¢ and s to obtain 


(3) 


= cosh k@ - sinh k@ 
2k 


and 


(4) 


] 
s(@) = = sinh k@ 


From x + ry =VN [c + rs) we find that x = WN c; y=Ws and 


(5) 1 k 
6=- sinh”? = 


k VN 


From (2) we also have [c(9) + rs(@)]" = c(n@) + rs(nO), a gener- 
alization of De Moivre’s theorem. For p=0, q=1 this reduces to 
[cosh 6 + rsinh 6]" = cosh n@ + r sinhn@ and we can obtain coshn@ 
or sinh n@ in terms of functions of cosh 6 and sinh @ by expanding 
and equating components. 


To obtain the relationship between c(@) and s() we substitute 
r, and r, in (2) and multiply the results. 


— 2 
= 4 (r, + + 


1 = ¢? + pes - 


To find the addition formulas for the functions c(6) and s(6) we 
substitute 6, and 6, in (2) and multiply. 


(rop/2)% - ,(r-p/2) 


[c(6,) + rs(6,)] [c(6,) + rs(@,)] = + + rslO, + 


= ¢(6,)c(6,) + rle(6,)s (6,) + c(6,)s (6,)] + (pr + q) s 


Equating components we obtain 


(6) c(6, + 6,) = c(6,) + qs(6,) s 


(7) + = s(9,)¢(6,) (6,) + ps@ ,) s (6,) 


" 
= 
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If in (6) and (7) we set 6, =- 6, =-@ we obtain 
l=c + qs(@)s (-6) 


0 =s (@)c + s (0) + ps(@)s (-0) 


from which 
-s(@) 
= 
+ pc) -ee%0) 


c (-6) = c(6) + ps(@). 
Using these results in (4) and (7) we have 
c(6, - 6,) = c(6,)c(6,) + pelO,) s (6,) - qs(O,) s 


s(6, -6,) = - s(6,)c(6,). 


Continuing, we are able to derive formulas analagous to the double 
and half-angle formlas. If we compute the modulus of the sum of two 
multiple numbers we obtain a generalized law of cosines. We are able 
to differentiate the functions c(@) and s(0) together with their 
inverses and finally we can show that the components of a function 
of a multiple number satisfy the partial differential equation 
Uyy - qUxx - pUxy = 0 which reduces to Laplace’s eq. ir tne trig. 
case and the wave eq. in the hyperbolic case. 

We illustrate one use of mltiple numbers by solving the Diophan- 
tine equation 

2? =u? + v?, 


let c= x+ry where r* =pr +q. If p= 0, q= -1, N(c) =x? + y? 
and =(x? ~ 3 xy?) +(3x4%y - y?)r. 


Using c = ye"? and N= y? eP® it is easy to show that [N(c)]” = 
N(c™), n = 1,2,3,... . Hence for n=3 we have 


(x? + y?)3 = (x4 Bxy?)2 + (3x2y y3)? 
from which 
= 
u 
v= 


This agrees with the general solution of elementary number theory. 


U.S. Naval Academy 


TEACHING OF MATHEMATICS 


Edited by 
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ect you, a interested, or questions 

e others to discu idli 


Angle of Inclination and Curvature 
David Gans 


The familiar notion of angle of inclination of a straight line, as 
presented in texts on analytic geometry and assumed in calculus books, 
seems so natural and perfectly straightforward that it comes as some- 
thing of a shock when we first learn that the notion involves a dis- 
continuity when applied to even the simplest of curves. If a straight 
line g meets the x-axis in a point P, the angle of inclination of g, 
we know, is the positive angle, less than 180°, whose initial ray 
extends from P along the x-axis in the positive direction and whose 
terminal ray extends along g in the direction in which y increases. 
The angle of inclination of a straight line perpendicular to the 
y-axis is taken either as 0° or 180°. 

To exhibit the discontinuity referred to, let us consider a very 
simple curve, say, one which is represented in some interval by a 
single-valued function y = f(x) possessing at least a first and 
second derivative. At a point where the tangent line to this curve 
is not horizontal it is clear that the angle of inclination of this 
line is defined, in accordance with the above definition, and, more- 
over, is continuous in the neighborhgod of the point. The angle of 
inclination of a horizontal tangent line at a point of inflection is 
not defined, the above definition being ambiguous for horizontal 
lines. Let us, then, agree to define it to be 0° at a point ofthis 
type on each side of which f(x) is an increasing function, and to be 
180° at a point of this type on each side of which f(x) is a de- 
creasing function. This agreement is reasonable since it achieves 
continuity for the angle of inclination at these points. 

There remain the horizontal tangents at points of relative max- 
imum or minimum. If we p-oceed along the curve and approach a point 
of relative maximum from the left, we find that the angle of inclination 
tends to 0°, whereas if we approach it from the right this angle 
tends to 180°. For a point of relative minimum the angle tends to 
180° or to 0° according as the point is approached from the left or 
the right. It follows that, no matter which of the values 0°, 18 is 
taken as the angle of inclination at a point of relative maximum or 
minimum, there will be discontinuity atthe point as far as the angle 
is concerned. 
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To sum up, using symbols, if @ represents the variable angle of 
inclination of the tangent to the given curve, then a is a function 
of x: @ = a(x), which is defined (or definable) for all values of x 
in the given interval, and is continuous for all these values except 
those corresponding to relative maxima and minima of f(x). 

The function a = a(x), of course, has no derivative at its points 
of discontinuity. A very large number of books on calculus, strangely 
enough, have overlooked this fact in their discussion of curvature 
in a plane. Differing somewhat in their definitions of curvature, 
these books all end up by computing it in the same way, namely, by 
finding the rate of change of the angle of inclination with respect 
to the length of arc. In the process they differentiate a= a(x), 
after writing it in the ambiguous forma = tan~'[f'(x)], then sub- 
stitute the result in the equat‘on da/ds = (da/dx)(dx/ds), and obtain 


(1) K = y"/(1 + y'2)3/2 


as the expression for curvature. They then proceed to use this formla 
to find the curvature at the vertex of the parabola 2% = day, at the 
end of the minor axis of the ellipse b2x? 4 a*y? = a2b?, at the point 
x = 7/2 of the curve y = sin x, and so forth, points for which their 
derivation of formla (1) is invalid because of thediscontinuity of a(x). 


Interestingly enough, the curvatures thus obtained at points of 
relative maximum or minimum are actually correct since formla (1) 
holds for such points despite the flaw in its derivation. This can 
be seen by noting that there is an alternative way of defining an 
angle of inclination which involves no ambiguity or discontinuity 
and which validates the faulty derivations without changing any of 
the formal steps. According to this alternative method, the angle 
of inclination, T, of a straight line with slope # is the principal 
value of tan~'s. Thus, the angle 7 is positive acute, negative 
acute, or zero according as@is positive, negative, or zero. Since 
T = tan ‘u is a continuous and differentiable function of u whenT is 
restricted to principal values, and f‘(x) is assumed to be dif- 
ferentiable throughout the given interval, so will T = tan~*[f‘(x)] 
be continuous and have a derivative, d7/dx = f"(x)/(1 +[f'(x)]?), 
throughout the interval. Substitution of this in d7/ds = (dr/dx)(dx/ds) 
gives formla (1). 

Fortunately, discussions of curvature can be improved without 
taking the somewhat drastic step of changing the definition of an 
angle of inclination. All that is necessary isthat the angle 7 of 
the alternative definition be used as a measure of the directim of a curve 
without calling it an angle of inclination. No special name for it seems 
necessary. Incidentally, this use of T would also represent a significant 
application of the principal values of the inverse trigonaetric functions, 
such applications being quite scarce in introductory courses in calculus. 


New York University. 


NOTES ON CIRCULAR AND HYPERBOLIC FUNCTIONS 


William S. McCulley 


1. One of the aims of instruction in calculus is to prepare stu- 
dents for solving diflerential equations. An important class of dif- 
ferential equations consists of ordinary equations, of second and 
higher order, with constant coefficients, which describe oscil lating 
and non-oscillating electrical and mechanical systems. The former 
are quite naturally described by oscillating functions, i.e., cir- 
cular sines and cosines. The latter are also quite naturally described 
by non-oscillating functions, i.e., hyperbolic sines and cosines. 

The usual syllabus for calculus courses places much more emphasis 
on developing and using the properties of the circular functions, 
relegating the hyperbolic functions to a coverage of two orthree les- 
sons. The closely parallel analytical prpperties of the two types 
of function render it feasible and instructive to treat them to- 
gether. Further, various relations between circular and hyperbolic 
functions are interesting to explore in themselves. 


2. As an aid to remembering many of the relations between circular 
functions the following diagram, seen originally in the preface to 
Wentworth & Smith’s Trigonometry, published in the early 1920's, 
serves as a useful condensation of a large amount of information. 


[\/\ 


Figure 1. 


It is easy to verify the following rules: 


2.1 a. Functions and co-functions appear to left and right, 
respectively. 
b. Derivatives of functions (on left) are positive. Der- 
ivatives of co-functions (on right) are negative. 
2.2 Reciprocal functions appear at opposite ends of diagonals. 
2.3 Each function equals the first function following divided by 
the second function following, in either direction. 
sin 6 
E.g-, tan 6 = 
cos @ 
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2.4 Each function equals the product of the two adjoining func- 
tions. E.g., cot 6=cos 6 esc 
2.5 The two functions at the top of the hexagon have finite am - 
plitudes; the other four become infinite. 
' 2.46 The Roman numerals indicate the pairing of functions in qu- 
adratic identities. 


Similarity of definitions permits a similar, but not identical, 


arrangement of hyperbolic functions, exemplifying all the above 
rules except 2.1 a. 


Figure 2. 


We shall see in section 5 how a comparison of Figure 1 and Figure 2 
yields further information about the graphs of circular and hyper- 
bolic functions. 


3. When we are given the value of one of the circular functions and 
wish to find the values of the other five functions without consult- 
ing tables, we may use the well-known device of reference triangles. 
A simple modification enables us to use a similar procedure for 
finding the values of hyperbolic functions (1). Suppose sinh 9= 3/5. 
Draw two right triangles. On the first of these triangles assign 
values to the altitude and hypotenuse according to the definition 
of the circular sine. On the second triangle, assign the same pair 
of values, but rotated clockwise one place with respect to sides. 
Thus we have 


(1) 5 (2) (34)? 


Figure 3. 


On the second triangle the length of the hypotenuse is obtained by 
means of the Pythagorean theorem. This value is assigned to the un- 
numbered side of the first triangle. The values of the remaining 
five hyperbolic functions may now be obtained from the numbers on 
the first triangle by applying, the definitions of the corresponding 
circular functions. 


4. The differentiation formulas for the circular and hyperbolic 
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functions and their inverses are assumed known. (See (2), (3), (4) 
“or details.) It is pertinent to observe that teaching and learning 
these formulas by pairing or “ stacking’’ them. E.g., 


cos sin 
D.(coshu) J sinhuDu, J =-1 (cire.); J = +1 (hyp.) 


tan 
D (are tanh u) > , upper function - upper sign, etc. 
+u 


5. Consider two parametrizations of the unit hyperbola x?-yel: 


x =cosh ?, y =sinh?, and 
(5.1) 


x=sec6, y=tan@., 


We thus establish a correspondence between the hyperbolic angle ¢ 
and the circular angle 6 expressed by 


(5.2) cosh = sec 9, sinh ¢ = tan 6. 


These relations can be symbolized by superimposing the hexagons of 
Figure 1 and Figure 2, from which alignment we can read off the re- 
maining four equivalences: 


(5.3) tanh sin 6, coth ¢ = esc 9, 


sech cos esch ¢ = cot 8. 


If we solve each of these, in particular (5.2), for 6 in 
terms of ¢, or for ¢ in terms of 6, we get 


(5.4) @ = arc tan (sinh ¢), ¢ = are sinh (tan 6), etc. 


It is clear from these relations that the domain (~©,©) for ? cor- 
responds to the domain (-',%m) for 9, and that ¢ = 0 corresponds 
to 6 = 0. We see that for corresponding values on the two domains 
the functions in{5.2)and(5.3)have equal function values, where they 
are defined, or else they both become infinite. Comparison of the 
evaluated derivatives of the functions in(5.2)and(5.3)shows that the 
corresponding functions have equal slopes at corresponding points 
of their domains. 

It can be verified by differentiation that d@/d? = sech d for 
each of the relations represented by(5.4a);differentiation of each 
of the six represented by 6.4b)leads to d?/d9 = sec 6. That dO/dp 
- dy/d@)~'may be seen from(5.2a, If we reverse these differentia - 
tions we obtain six equivalent integration formulas for sech ? and 
and for sec 

If we wish to derive such formulas by formal integration, we may 
may proceed as fol lows: 
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d¢ cosh? dp 
(5.5) sechddd JS; sinh’? + arc .tan (sinh $) +C, 
etc. 

It is interesting to compare the number of operations necessary to 
evaluate antiderivatives of the form(5.4b)and(5.5)with the number re - 
quired for evaluating the standard antiderivatives. In the case of 
(5.4b)the standard formula is {sec Od0 = In(sec @ + tan 6)+C, 
which requires three references to tables and one addition; (5. 4b) and 
(5.5) can be evaluated with two references, or, on a vector slide 
rule, by one setting of the indicator. 

In his work published in 1830 Gudermann used the relation(5.5) 


to define the function he named the longitude. Cayley later renamed 
it the gudermannian. Thus 


(5.6) gd od =f sech = arc tan (sinh C, 
and its inverse 
(5.7) =f sec Od 6 = arc sinh (tan 6) +C. 
From the discussion following(5.4)we find gd(0)=0, gd(tm)=+ Kn, 


6. From the correspondence relations cos 9 = sech ¢, sin 0 = 
tanh d, we may parametrize the unit circle x* + y? = 1 in terms of 
the hyperbolic angle ¢ by x =sech¢, y=tanh¢. Then the area of 
the portion of this circle in the first quadrant is given by the 
line integral 


(6.1) x y sech d tanh ¢ 
A= % f | ip = % f dp 
dx/dp dy/dd -sech¢tanh¢ sech* ¢ 


ech = % ed = 0) = 
Similarly, the length of arc is given by 


6.2) s =f" tanh’ + sech*¢) tap sech ¢ap =gd | - Xn, 
0 0 


7. The functional correspondence established in (5.2)and(5.3) gives 
us another way of integrating sec’9, hence odd powers of séc 0, as 


follows: 
sec = ff sech = ff 
= % (sinh cosh + + C 


(7.1) 


= % (tan @ sec 6 + arc sihnh (tan 6)) + C. 


The same functional correspondences may also be used to establish 
hyperbolic analogs of the Wallis formulas, thus 
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tn © © 
f sin’? do = J = sini? cosh" dd; 
(7.2)° 


cos"@ dO = ff sech™' dp =f. dd, etc. 
0 0 0 


8. As an example of the compact expression obtainable by combining 
circular and hyperbolic functions, consider the differential equation 


(8.1) Yee + 9%, + Oy = 


where a and b denote constants. The auxiliary equation is 


of which the roots are 
= 


The particular solution is given by 


(8.1) i 
y= e**(Acosh pt + Bsinh pt), 
2_& 
where a=-%a, 8=% 4b)]*, with the 
rule: 


If J=+1l, use hyperbolic functions; 1f J=-l, use circular functions. 


9. Additional analogies and generalizations involving circular and 
hyperbolic functions include the fol lowing: 


a. The hyperbolic analog of de Moivre’s theorem 
(cosh + j sinh = cosh n¢+J sinh nd, where j* = +1; 
b. The generalization of circular and hyperbolic functions to 


oscillating and non-oscillating Bessel functions respective! y. 


c. The generalization of circular and hyperbolic functions to 
the elliptic functions. In particular, the Jacobian elliptic 
functions sn(u|m), cn(u|m) satisfy, for limiting values of 
the parameter m, the following equations (see [5] ) 


sn(u|0) = sinu, sn(u|1) = tanh u, 


en(u|0) = cos u, en(u|1) = sech u. 
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The Derivatives of the Trigonometric Functions 


M. J. Pascual 


The following method of arriving at the derivatives of the trig- 

onometric functions might appeal to some, even if for the sole 

reason that it differs from the conventional method employing the 
sinh 

- A novel feature of this derivation is the fact that the 


derivative of the sin @ is arrived at last and then as a consequence 
we may obtain the above limit. The derivation will be for © lying in 
the first quadrant and with AO positive. Obvious modifications can 
be made for other quadrants and A@ negative. 


R(x,y + Ay) 


We start with the tan @ and cot @. By definition 


6 + AD) tan? 


y 


6 


where x 7 0, 0677/2. Similarly we obtain 


where y # 0, @ #0. Now to obtain the lim x A from the 
above diagram it is easy to see that ti 


area OPS < area OPR < area OAR 
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Ss 
P(x,y) 
x 
x x ] 
A 
8) lim x = 
646 
|_| 
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% (x2 + y?) < dy < % [x? + (y + Ay)?] AO 
hence 
x = x74 y?, 
There fore 


x? 4 y? 


D (tan 6) = ——— for x#0, 64 n/2, 


x2 
D (tan 8) = sec? 6 


+ y? 
D (cot = ———— 


y #0, 670 


D (cot 8) = = 


By using the identities relating the trigonometric functions we 
easily get 


(1) D (sec 9) = sec 9 tan 9 and D (cos 9) = - sin 9 for 64n/2 
(2) D (esc 6)=-csc cot and OD (sin 8) = cos 9 for 


To verify that the formlas for the derivatives of the sin@ and the 
cos 6 hold for these exceptional values of 6 as well, we have by the 
quotient formula that 


6 i - sin 0 6 
cos D(sin 9) - sin 8 D(cos @) tee OF 


cos?@ 


cos D(sin + @ 


sec? 6 = 
cos*6 
solving for D(sin 6) we get D(sin 6) = cos @ for 6 7# n/2. Combin- 
ing this with (2) we may say that for 0 < 6 < 1/2 


D(sin 8) = cos 8 


Similarly by using D(cot g@) for @ 7 0 we obtain D(cos 9) = - sin @ 
for €# 0 which combined with (1) yields 


D(cos 8) = -sin for 0< 6< n/2. 
sin(6th)-sin 9 


Finally to evaluate the lim = we have D (sin 9) = Lim 


sinh 


which for 6 = 0 gives us_ cos 0 = lim 
h-70 A 


Siena College, New York 
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ON A CHARACTERIZATION OF ORTHOGONALITY 
Waleed A. Al-Salam 


D. Dickinson [1] gave the following theorem: 


Theorem: For a set of polynomials P, (=), n= 0,1,2,... where each 
P, 1s of degree precisely n, to have the property 


a, = Pa. pyl-x) + ,(-x)p, (x) 


for some a, ast necessary and sufficient that 
there exists a recurrence relation of the form 


p,(x) = xB, Pn + Cop, s(x), Cc, #0 (n=2,3 
and that p,(0) # 0. 


He also proved the following characterization of orthogonality 


Corollary: For a set of polynomials p, («), n = 0,1,2,... where each 
Pp, 18 of degree n, to be orthogonal, it is sufficient that there 
exists a relation of the form 


= Paw + Pal -x) 


where (-])" a,a, < 0, fer 


In this note we obtain a more general criterion of orthogonality 
which actually reduces to Dickinson's criterion as a special case 
and further includes the classical polynomials. We also obtain the 
orthogonality of a class of polynomials associated with the ortho- 
gonal polynomials. We first prove the following theorem: 


Theorem 1. Given two sequences of polynomials { f.. (x)} and 
{g (x)} n= 0,1,2,..., where jf. {x) and yg, (x) are of degree n and 
n-1 respectively, fo # 0, 7 9, go = 0, then the necessary 
and sufficient condition that there exists a relation 


is that {f,} and {y,} satisfy a recurrence relation of type 


(2) (A, x + B,) u(x) + Cu 


(x), (C40, n=1,2,...) 


Proof: If {f (x) } and {g,,(x)} satisfy (2) then it is easily seen 
that 


4] 


——— 
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= (-1)"C,C. 
7 0 for all n = 1,2,... 


Now assume that (1) holds. Then by Casorati’s theorem {f,} and 
{g,} are independent solutions of the difference equation 


(x) 


Unel 


u, (x) f,(*) 6, (x) 
(x) 


(x) 


n+1 Bnei 


By expanding we find 
x Sas 18n-1 


a, 
u,(x) - —— u(x). 


n- i 


1 (*) = 


Now consider F,(x) = f,.s8n41 If we eliminate 
by means of (1) we get 
fn 
fae 1%n 


fn 


But a, does not involve x, and f,,, is of degree n + 1, then we get 
that P must be of the first degree. Hence 


(x) 


a 
(x) = u(x) -——u 


Une nel 


This completes the proof of the theorem. 
In particular if we have g(x) = f,(x) + (-1)"**f,(-x) then the 
expression in the L.H.S. of (1) reduces to 


n+1 


and similarly the expression F, reduces to 


| | 
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from which we see that F.(0) = 0 and hence B= 0. This shows that 
given a sequence {f (x)} f, isa of degree n,f 9 #0, 
then the necessary ‘and sufficient condition that there is a relation 


f=) + = a, 0 
is that {f ,) satisfy a difference equation of type 


(3) = (x) + Cu (x), 0 


which is Dickinson’s theorem. Another way of stating this is that 
the necessary and sufficient consition for {f, (x)} and {(- 1)"f (-x) } 
to be two independent solutions of (3) is that there is a relation 


for some a. 


Now the following characterization of orthogonality can be proved 
by using Favard’s theorem [2) and the previous theorems. 


Theorem 2: Given a sequence {f,(x)} where f, is a polynomial of 
degree n and fy 7 0 then the necessary and sufficient condition 
that {f, (x) } be orthogonal is that there exists another sequence 
of polynomials {g,(x)} where g, is a polynomial of degree n - 1, 
Bo = 9, AO such that 


(2) = a, 40 and aoa 0. 


It is well known that if {f (x)} satisfy the recurrence relation 
(2) then we have for each p = 1,2,3,... 


fae p(®) = (x) + 


where T(")(x) and re are polynomials in x of degrees p wd p-! 
respectively and where (z) = I, T{")(x) +B, = 0, 
and (x) = C, 7 0. Let us call T(")(x) and for each fixed 


n the Loinme | polynomials of the and kinds associa’ 
with f, (x) respectively. Then we can prove the following 


Theorem 3: Given a sequence of orthogonal polynomials {f,, (x) }where 


deg f, = n, and fy A 0, then the Lommel polynomials associated with 
each Tixed n are also exrtiigaint. 


Proof: Since the {f (x)} are orthogonal, they satisfy a difference 
equation of the second order. 


n 
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(x) = u(x) - Cy, 


Unet 
and there exists another sequence {g,(x)} such that 
fy (*) fag) = #0 where a,a,> 0. 


Now consider the Lommel polynomials associated with f, (x). We have 
after Palama [4]. 


(4) a, (x) = fy Bay p(*) Byes 
and 
(5) a, 50") (x) = fay (*) - f,(*) 


Hence the expression 
a, h(x) Mx) - Thm (x) 


after substituting (4) and (5), reduces to a,_ ,a,,, and by theorem 
2. thie fomen is established, i.e., the two sequences (ro (x)} 
and ")(x)} p = 0,1,2,3,... are orthogonal. 

The " Cathenien expressions of type (1) are known for the classical 
polynomials ( see Toscano, [3] ) 


+n + 1) 
FOO) ) (x) (x) = 


(n+ 1)! (a+ 1) 
(2a), 
(n+1)! 


(x) PCN x) (x) = 
H (x) - H, G(x) =n! 


where L (Xx), P(2)(x), and (x) are the Laguerre, the ultrasherical, 
and the Hermite polynomials respectively, whereas FO) (x), RO) (x), 


G(x) are the polynomials associated with them. 


eee 
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ON NATURAL BOUNDARIES OF A GENERALIZED 
LAMBERT SERIES 


Francis Regan and Charles Rust 


1. Introduction. In 1932, Feld [1]! introduced the series 
@ 5 n 
Fiz) = 5 (1) 
l-a 2" 


and showed that any analytic function can be represented in this form 
plus a constant. Doyle [2] determined the regions of convergence, ex- 
pansion in power series and the inversion of power series into a 
generalized Feld series, which included (1) as a special case. 

It is the purpose of this paper to investigate the Feld series for 
natural boundaries. In establishing natural boundaries for this series, 
conditions on the sequences {a} and {b, } of (1) are determined. Re fore 
obtaining these conditions three theorems dealing with natura! bound- 
aries of the Lambert series 

n 
L(z)= = (2) 


mel) 2" 


are presented and are used in broadening the conditions for natural 
boundaries for (1). 


2. As an immediate consequence of the work of Knopp [3], three theorems 
dealing with natural boundaries of the Lambert series follow. 


Theorem 1. If {c,} is a null sequence of positive real numbers, 
then the unit circle is a natural boundary of the 
Lambert series (2). 
Since {c,} is a null sequence, then 


for a denumerable set of k’s. Hence the unit circle is a natural 
boundary of (2). 


We come next to 


1. The numbers appearing in [] refer to reference given in bibliography 
at end of paper. 


45 


n 
= Cc, 
1 
lim = 0, 
n 
kr<n kr 
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Theorem 2. If {c,} is a real sequence such that 0<c <c¢,< C<®, 
then the unit circle ts a natural boundary of the function represented 


by the Lambert series (2). 


Using Knopp’s method, we let z, = io with (k’,k) = 1, be any 


rational point on the unit circle and then from the Lambert series, we 
have 


(2) 


Multiply both sides of (3) by (1 -2z/z,) and take the limit as z ap- 


proaches z) radially. We examine the two limits on the right. Hence 
k 
Lim (1 > 
vj] 


which, if C* = y, we see that : 
v 


Chey? 
(1 -2z/z9) = 


On the other hand 


leall2 


- 


| (] - 


When n is not a multiple of k, |] - z"|> h>0, where h=1 when k=2 
and h = sin 27/k when k > 2. Hence 


(1 - 2/2.) 
_ 2” h h 


Since (4) is unbounded and (5) bounded, we have 


n 
c,2 


lim ja = 
1 - 2* 


46 
n kv n 
@ Cy Z 
——= —., = 
nzi ] vel -z* ken 1 
lim 1 
c y ] 
(5) 
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Hence the unit circle is a natural boundary of the function represented 
(2). Lastly, we have 


Theorem 3. If {c, } ts a real sequence such that 0<c¢,<¢ C <@ 
ond if for a denumerable set of positive integers k, the series 


E ey, /kv is divergent, then the unit circle is natural boundary of 
v= 


the function represented by the series (2). 


As before we consider the two limits in (4) and (5). From (4), we 
easily get 


lim 


((4), p.177), and it follows easily that the second is bounded. When- 
ever, this is true for a denumerable set of k’s,then the theorem follows. 


3. Before discussing the general Feld series, we shall state some 
theorems which are special cases, the proofs of which are obvious. 


Theorem 4. If a = a” and b= 1, then the Feld series represents 
a function which has'the circle |z| = 1/|a| as a natural boundary. 
The analog to Knopp’s extension of the Frane] theorem is 


Theorem 5. If the radius of convergence of = a"b z” is greater 

than or equal to 1/ |a|, and if the numbers 6b, are such that for a 

certain k, the series p3 with 1 7 0, 1, ..., k-1, converge if 
kv +l 

for such a ka relatively prime integer k', we set = 2" 


then for radial approach 
2" 


lim - 2/2.) 


1 a™z" 


If for a denumerably many k = (b,,/kv) 7 0, then the circle |z| = 1/|a| 
is a natural boundary of the function represented by (1). 


Theorem 6. If the radius of convergence of = a®b z” is greater 
than or oqust to 1/|a|, and the b, are real and positive, the series 
a"b_2z” 
n 


= ———— _ represents a function which has the circle |z| = 1/\a| as 
1 - a"z" 


a natural boundary when any one of the following conditions on the b, 
is satisfied 


v0, 

o “ky m Lin — y 

vel py 
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(a) when the numbers b, forma null sequence, or 


(b) when the numbers b, satisfy the inequalities 
0<bsb SB<®, or 


(c) when the numbers by, satisfy0 <b SB<o@ 
and the series 5 
(b,,/kv), 


diverge for a denumberable set of positive integers k. 


Now let us consider the series when a, =a and 6, = 1. Only when 
ja| < 1 will be investigated for if a=] the Feld series (1) reduces 
to the Lambert series and if |a| > 1, by means of the transformation 

len n 
w 
az = w, we obtain where <1, when n> 
he 
az" 
Lemma. If |a| <1, then = => for 1. 


1 - az” ic 


This lemma is easily established by using Weierstrass theorem on double 


series. 
az” 


Theoren 7. If la| < 1, then the series = ———, represents a 
1 -az” 
function having the unit circle as a natural boundary. 


From the lemma above, we readily obtain 


a™z" 


az" 
lim { (1-2/z,.) 5 lim 


9 1 n 


where z, is any rational point on the unit circle. Hence the unit 
circle is a natural boundary. 


Returning to the general Feld series (1), we impose two conditions 
on the numbers a,. Firstly, | < 1 and secondly, Lim | 


Now assuming that the radius of convergence of = a"b z"is greater 


than or equal to unity, the Feld series is uniformly convergent for 


a,b,2” 


|z| < 1. Expending in a power series, we have ———— = A,2", 
1 - 
n/d 
where A, a) 


the sequence {a,b} may be expressed in terms of the sequence {A,}. This 


b,- By means of a generalized Moebius function {2], 


function is defined as follows. 
S(Q1) =1 


GENERALIZED LAMBERT SERIES 


(u/d)-1 
d/u S(d) =f foru > 1. 


We can show that within the region of convergence of > a,6,z" and for 


a 


<1, 


=ZA,z", where a,b, = as folloss. 
n 


d/d 
Multiply both sides of A, - E oy, ‘ba by S(n/d) and summing over 
1 


all disvisors of n, we obtain 


n/(d,d) 
>» ‘d =z = 1 
bin S(n/ YA, a, ba, d, 


n/(d,d)-1 
d 1 
“dy “dy dy aq, 


2 S(d) = 
n d (a 
1 
n/(d 


=< b, = S(d)a 
1 


«1 
But let n/d = u, the coefficients of each a, become 
u 


which from the definition equals zero for all u = 1 and equals one 
when u = 1. But when u=l, n/d, = 1 ord, =n. Therefore 


c 
Now, expand the Lambert series (2) in a power series, we get = —a*_= 

A,z", where A, = . 


If we now let this A, be the coefficients of the power series we ob- 


a,b,z" 
tained by expending the Feld series, we will have 5 ————— = 5 
1 -a,2" 1 - 


where a,b, = We now conclude that if |c,|<C <@ 
and Ja, | < 1, then the Feld series and the derived Lambert series re- 


present the same function within the unit circle. Since these two series 
represent the same function within the unit circle, if we multiply both 
by the same factor and take the limit as a point is approached radially 
om the unit circle, we will obtain the same limit. Hence using Knopp’s 
theorem for the Lambert series, we easily eatablish 
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Theorem 8. If the numbers a and b, of the Feld series are such 
that la, | < and ab, = C4, , where for an integer k, 


if for such a k and a relatively prime k', we set Zy= e : , then for 


radial approach 


Ch. 


] «46 2" 
n 


c 
v 
If = {> 79 for a denumerbale set of integral values of k, the unit 


circle will be a natural boundary for the function represented by the 
Feld series. 


If we make use of Theorem 2, we conclude 


Theorem 9. If the numbers a, and 6, of the Feld series are such 
that |a,| <1 and = S(n/d) 4, where SC <a, 


then the unit circle is a natural boundary for the function represented 
by the Feld series. 


From Theorem 3, we have 
Theorem 10. If the numbers a, and b, of the Feld series are such 
that < 1 and ab, = where 0< c, <C < wand 


c 
k 
= _ is divergent for a denumerable set of k’s, then the unit dircle 


is a natural boundary for the function represented by the Feld series. 


BIBLIOGRAPHY 


1. Feld, J.M., The Expansion of Analytic Functions in a Generalized Lambert 
Series, Annals of Mathematics, Vol. 33 ( 1932), pp. 139-43 

2. Doyle, William, A Generalized Lambert Series and its Moebius Function, 
Annals of Mathematics, Vol 40 (1939), pp. 353-59. 

3. Knopp, Konrad, Ueber farterteche Reihen, Journal fur Mathematik, Vol. 142 

(1913), pp. 283-3 
4. -——, —, Theory and ieee of Infinite Series, London, Blackie and 
Sons, 1946. 


St. Louis University, St. Louis, Mo. 
Xavier University, Cincinnati, Ohio. 


PROBLEMS AND QUESTIONS 


Edited by 
Robert E. Horton, Los Angeles City College 


Readers of this department are invited to submit for solution problems be- 
lieved to be new and ochdons matter questions that may qtico in study, in 
research, or in extra-a ¢ situations. Proposals should ag ee by 
solutions, when available, and by any information that aril a e editor. 

inarily, problems in well-known textbooks should not be submitted 

Solutions should be submitted on separate, signed sheets. Figures should 
be drawn in India ink and twice the size desired for reproduction. 

Send all commnications for this department to Robert £. Horton,Los 
Angeles City College, 855 North Vermont Ave.,Los Angeles 29, California 


PROPOSALS 


313. Proposed by Sidney Kravitz, East Paterson, New Jersey. 


Show that Christmas falls on a Sunday more often than once every 
seven years. 


314. Proposed by J.M. Gandhi, Lingray College, Belgaum, India. 
Show that = 2(.C, - ,C,) + - + - 
the last term being n(\C_) when n is even and -(n-1),C whennis odd. 


315. Proposed by P.D. Thomas, Eglin Air Force Base, Florida. 
Under the restriction f(c - x)= f(x) show that 

Sjla + bx) f (x) dx may be written (a + (be)/2) [£f(x)dx. 

316. Proposed by A.K. Rajagopal, Lingraj College, Belgaum, India. 


Prove that (1 + cos n9), n an integer, has a factor (1+ cos 6) 
if and only ifn is odd. 


317. Proposed by Ben K. Gold, Los Angeles City Colleye. 

Prove that (e + > (e x)f¢**) for x <e, 
318. Proposed by Chih-yi Wang, University of Minnesota. 

log x 

Evaluate lis 
319. Proposed by Barney Bissinger, Lebanon Valley College, Pennsylvania 
Cos Cos 26 Cos(n-1)¢ 
+ 


sin nO 


for those values of @ for which the terms are defined. 


51 


MATHEMATICS MAGAZINE 


SOLUTIONS 
Late Solutions 
274, 286, 288, 289, 290; B. Keshava, R. Pai, Mangalore, India 


286, 287, 288; Gene B. Parrish, Durham, North Carolina 
291; A.K. Rajagopal, Lingray College, Belgaum, India 


A Planar area 


292. (January 1957) Proposed by Eugenio Calabi and Chih-yi Wang, 
University of Minnesota. 


Find the area of the region in the real xy plane such that 
|sinh x sinh y| <1. 


I. Solution by Gene B. Parrish, Burham, North Carolina. From the 
symmetry of the area for which -1< |sinh x sinh y|<1, the area in- 
tegral may be confined to the first quadrant, with x ranging from 0 
to @ and y ranging from 0 to sinh~‘csch x = In[(e* + 1)/(e* - 1)). 
Denoting by A the first-quadrant portion of the area, the total area 
sought is 


@ 
4A=4 Jf, In [(e* + 1)/(e* - 1)) de 
Letting (e* - 1/(e* + 1) =u, then 
w du/(1 u?) = 1? 


where the last integral is given in Table 108 of Bierens de Hans’, 
Nouvelles Tables D’Integrales Definies and in Edwards’, Treatise on 
Integral Calculus, Vol. II, p. 249. 


II. Solution by Chih-yi Wang. By symmetry, it suffices to find the 
area in the first quadrant. Let A be the total area of the required 
region. By using the relation sinh=‘u = log (u + /1+u2 ), we have 


A/4 = J, sink (esch x) dx = J, log (csch x+coth x) dx 


@ 1 + cosh x aw+il 
= l dx = 
0 sinh x e* -1 
(Pierce 520) 
Therefore, A = 77? 
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III. Solution by Eugenio Calabi. Let x = sinh™'(tanu cos v), 
y = sinh~'(cos u tan v). Since the Jacobian 


(x,y) 
(u,v) 


J 


iF 


the area in question is the same as that of |sinusinv | <1, 
but we know this relation holds for |u| < 7/2, |v| < 7/2. 
Hence the answer is 77”, 


An Euclidean Asyaptotic Construction of e 


293. (January 1957] Proposed by Raphael T.Coffman, Richland, Washington. 


Given a line of unit length, construct geometrically a line of 
length (1 + 1/n)”, where mis an integer. 


Solution by Leon Bankoff, Los Angeles, California. Consider the 
more general problem of constructing geometrically a line of length 
(1 + a)", where n is an integer and a is the length of a constructible 
line segment. 

Extend the given segment AB > f so that BC = a and AC = 1 + a. 
Erect perpendiculars BB’ end CC’ to AC at B andC respectively. 
Describe arc (4,AC) cutting BB’ at D, and extend AD to cut ’ at &. 
Then AE = AD + AC/AB = AC*. Now describe arc (A,AE) cutting BB’ in 

_F, and extend AF to meet OC’ at G. Then AG = AF - AC/AB = AE + AC =AC?. 
Repetition of this procedure leads in an obvious way to the construc- 
tion of (1 + a)" for any desired value of n. For the special case, 
(1 + lfm)", let a = I/n. 


Ase = lim (1 + 1/n) we have an asymptotic construction of e. 


Also solved by Dermott A. Breaalt, Carnegie Institute of Tech- 
nology; David J. Cartmell, The College of Wooster, Wooster, Ohio; 
Howard Eves, University of Maine; B. Keshava R. Pai, Mangalore, India; 
Chih-yi Wang, University of Minnesota and the proposer. 


Concyclic Nine Point Centers 
294. (January 1957) Proposed by N.A. Court, University of Oklahoma. 


The nine-point centers of the four triangles formed by four con- 
cyclic points taken three at a time lie on a circle. 


I. Solution by Sister M. Stephanie, Georgian Court College, New 
Jersey. Let the circle on which the four points lie be a unit circle, 
and the four points be, in a system of complex coordinates, t,, t,, 
t, and t,. The nine-point circle of triangle t,tpt, will then be 
z= s,/2 + t/2 where s,; = t, + t, + t,. Its center is at s,/2 and 
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its radius in one-half the radius of the unit circle. The nine- 
point circles of the other triangles are found similarly. 

Consider the circle z = Sy/2 + t/2 where Sy = ty + tp + t3+ty- 
This circle is the locus of the centers of the four nine-point 
circles, for the substitution of a properly chosen unit vector, 
e-g-, -ty yields each of the centers in turn. It has center at 
S,/2 and radius equal to one-half the unit circle. 

Note: The point S /2 is the point of intersection of the four 

nine-point circles since it lies on each; it is also the center of 
the equilateral hyperbola which can be drawn through the four 
given points. 
II. Solution by B. Keshava R. Pai, Mangalore, India. The radius of 
the nine-point circle of a triangle is equal to half of the cir- 
cumradius of the triangle. Therefore, the radii of the four nine- 
point circles of the four triangles are equal, since the four points 
are concyclic and so the circumcircle, is one andthe same, viz 
the circle passing through the four points. 

Now, the centre locus of a family of Rectangular hyperbolas 
passing through the vertices of a triangle is the nine-point circle 
of the triangle. Consider each of the four triangles separately. The 
four nine-point.circles are the four centre loci of the four sets 
of Rectangular hyperbolas. Now, through any four points there passes 
a Rectangular hyperbola. The Rectangular hyperbola passing through 
the four given points belongs to all the four sets and hence its 
centre lies on all the four loci. Hence, the four loci, i.e., the 
four nine-point circles pass through a point. 

Since the four nine-point circles have common radii and pass 
through a common point, the centres of them lie on a circle with 
the common point as the centre and half the circumradius as radius. 


Also solved by J.W. Clawson, Collegeville, Pennsylvania; Huseyin 
Demir, Kandilli, Eregli, Kdz, Turkey; Howard Eves, University of 
Maine; A.K. Rajagopal, Lingraj College, Belgaum, India; Chih-yiWang, 
University of Minnesota and the proposer. 


Radix Identities 


296. January 1957 Praposed by P.A. Piza, San Juan, Puerto Rico. 
Prove that the following equalities 
385 + 439 + 547 = 3467 + 475 + 529 
385? + 4392 + 5472 = 3672 + 475? + 529? 
are true not only when the six distinct 3-digit numbers are considered 
to belong to the decimal system of numeration, but also when they 


are regarded as belonging to any system or scale of numerical 
notation with any base greater than ten. 
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Solution by James A. Painter, I.B.M. Corporation, Endicott, New 
York. A three digit number abc in any base, say r, represents 
ar? + br +c. Hence the problem is to show: 


(3r?2 + Br + 5) + (4r 24 3r +9) + + + 7) = 
(3r2 + 6r + 7) + (4774+ 7r +5) + + 2r + 9) 


(1) 


(3r? +8r+ 5)?2+ (4r? + 3r + + + 4r + 7)? = 
(2) 
¢3r2 + 6r +7)? + (477+ Ir + + (Sr? + 2r + 9)? 


for r 2 10. 


Collecting like powers of r gives 3 equations from equation (1) 
and 5 from equation (2). 


§+9+7=7+5+9 

(8 +3 +4)r = 2)r 

(3 +44 5)r? = (3 +4 + 5)r? 

25+ 81 +49 = 494 254 81 

(56 + 54 +80)r = (84+ 70 + 36)r 

(64 +3049 72416 +70? = (36 +49+4+4 42+40+90) r? 
(40 + 24+ 48)r2 = (36 + 56 + 20)r? 


(25+ 16 +9)r*= (9 +16 + 25)r* 
for r 2 10. 


Obviously this set of equations hold for any r. The base of a 
positional notation system must be greater than the value of any 
symbol used in the system. Since 9 is used in equation (1) and (2), 
r> 9. That is r 2 10. 


Also solved by Haward Eves, University of Maine; Harry M. Gehman, 
University of Buffalo; Erich Michalup, Caracas, Venezuela; Wahin 
Ng, San Francisco, California; B. Keshava R. Pai, Mangalore, 
India; Chih-yi Wang, University of Minnesota and the proposer. 


Construction of an Equilateral Triangle 


297. [January 1957] Proposed by Dewey Duncan, East Los Angeles Junior 
College. 
Given a point, a line and a circle in a plane. Construct an eq- 
uilateral triangle having a vertex on each of them. Determine the 
criterion for the existence of such a triangle. 


Solution by Howard Eves, University of Maine. Designate the given 
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point, line, and circle by P, L, and C. Let L’ be obtained by ro- 
tating L about P through an angle of +40°. Suppose L’ cuts C ina 
point Q. Then the perpendicular bisector of PQ cuts L in R such that 
triangle PQR is equilateral. 

A necessary and sufficient condition for a solution is that L’ 
intersects C. Because the angle of rotation may be +40°, there may 
be as many as 4 solutions to the problem. It is easy to construct 
examples showing exactly 3, 2, 1, and no solutions. 


Also solved by J.W. Clawson, Collegeville, Pennsylvania; Huseyin 


Demir, Kandilli, Eregli, Kdz, Turkey; H.M. Gandhi, Lingraj College, 
Belgaum, India and the proposer. 


an Invariant Curve 


298. [January 1957] Proposed by Huseyin Demir, Kandilli, Bolgesi, 
Turkey. 


Let y = f(x) be a curve with the following properties 


a) f(x) = fl-x) 
b) f(x) >0 for x>0 
c) fx) = 0 


Determine the weight per unit length w(x) at the point (x,y) such 
that when the curve is suspended under gravity by any two points 
on it, the curve will keep its original shape. 


Solution by K.L. Cappel, Philadelphia, Pennsylvania. Assume the 
curve to be suspended at two arbitrary points A and B. Let the 
weight between A and the y intercept of the curve be W. Then at 
A, the tension in the curve can be resolved into vertical and hor- 
izontal components so that W/H = tan Gor W=H dy/dzx. 

Now assume the right point of support to be moved fromA to 
A‘. If the curve is to retain its shape, there must be no change 
in the forces at A. This can only be the case if H is a constant. If 
ds is the length of the segment AA’, and dW is its weight, then the 
weight per unit length will be 


dw d °y/dx? 


| dy.2 7 


which can be satisfied by any curve obeying the given conditions. 
This problem is analagous to the problem of finding the optimum 
shape of a masonry arch, when the material of the arch is the only 
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load to be supported, and it is desired to have the thrust load act 
along the neutral axis in order to eliminate bending moments. 


Also solved by the proposer 


QUICKIES 


From time to time this department will publish problems which may 
be solved by laborious methods, but which with the proper insight may 
be disposed of with dispatch. Readers are eegee to submit their 


favorite problems of this type, together with the elegant solution 
and the source, if known. 


Q 204. If A, B, C and D are vectors and [ABC)is a scalar triple 
product, prove that 

[BCD] A - {CDA] B + [DAB] C - [ABC] D = 0 
[Submitted by M.S. Klamkin.] 


Q@ 205. An airplane with an airspeed of 100 mph, flying into the 
wind, passes over a point just as a lighter than air balloon is re- 
leased. After a while the plane turns down-wind and overtakes the 
balloon 8 miles from the point of release and one half hour after 
it was released. Assuming constant wind and no time lost in turning, 
how far did the plane fly before turning? [Submitted by Richard K. Guy.] 


Q 206. If x, y and z are positive and if x + y + z = 1 prove that 
-1)Q)/y - 101/z-1) 28 [Submitted by M.S. Klamkin) 


q 207. Find the Arithmetic Progression for which the nth term+ ath 
= (n + m)th term. [ Submitted by B. Keshava R. Pai.) 


208. If f(x) = f(x +1) = f(x +72) and f(0) find f(x). 
[ Submitted by M.S. Klamkin ] 


Q@ 209. Find the area between y = x? and y*= 32 «x. 
[Submitted by John M. Howell} 


Q 210. If f(x) can be integrated in finite form, show that the in- 


verse function f~'(x) can also be integrated in finite form. 
[ Submitted by M.S. Klamkin] 


ANSWERS 
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From the Preface to: 
“THE TREE OF MATHEMATICS” 


There is a great deal being published about mathematics these 
days, and that is fine; but this book is mathematics in the sense 
that it presents the epitomes of the main branches of the subject 
beginning with high school algebra and extending far into graduate 
work. 

When writing this treatise the authors gave great attention to 
making it both meaningful and understandable, an art in which most 
of them are pastmasters.The practice was to start on ground familiar 
to everyone and construct a highway, free from road blocks, through 
the wonderful world of mathematics. 

This book is, for the most part, a response to requests for 
source material from two classes of people: those who need an ever 
increasing knowledge of mathematics in their jobs, such as en- 
gineers; and others who have gone little if any beyond arithmetic 
and either need more mathematics or just want to know “what it is 
all about.” Inherent in a satisfactory response to these requests 
is an answer to the needs of teachers and students of mathematics 
and related subjects who desire to extend their horizons by home- 
study, and to do so in a minimum of time. 

In terms of class-work, the first seven chapters and perhaps 
selections from later ones would constitute about a one semester : 
survey course in high school, while by passing swiftly over the 
first three or four chapters the entire volume could be covered 
in a two semester survey course in college. 

But classwork is always slower than effective individual study. 
The latter could reduce the above periods to months or possibly 
weeks by study at home during the evenings. 

Due to its broad coverage, this treatise offers an excellent 
opportunity for “reading up” on specific topics. By looking a topic 
up in the index you will find references to the page or pages where 
it is discussed or its meaning depicted by usuage. Of course,a 
similar procedure is possible with textbooks, but it would require 
some two dozen of them to cover all the topics in this treatise. 

The Tree Of Mathematics, containing about 350 pages with 85 
cuts and pleasing format and typography, will sell for the moderate 
price of $5.50 if cash is enclosed with order or $6.00 if billing is 
required. This book will be ready for delivery in about 60 days. 
Orders will be filled in the sequence of their arrival. Address: 


The Digest Press, Pacoima, California 
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